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Abstract — The exponential synchronization rate is addressed 
for Kuramoto oscillators in the presence of a pacemaker. When 
natural frequencies are identical, we prove that synchronization 
can be ensured even when the phases are not constrained 
in an open half-circle, which improves the existing results in 
the literature. We derive a lower bound on the exponential 
synchronization rate, which is proven to be an increasing function 
of pacemaker strength, but may be an increasing or decreasing 
function of local coupling strength. A similar conclusion is 
obtained for phase locking when the natural frequencies are 
non-identical. An approach to trapping phase differences in an 
arbitrary interval is also given, which ensures synchronization 
in the sense that synchronization error can be reduced to an 
arbitrary level. 

Index Terms — Exponential synchronization rate, Kuramoto 
model, pacemaker, oscillator networks 



I. Introduction 

The Kuramoto model was first proposed in 1975 to model 
the synchronization of chemical oscillators sinusoidally cou- 
pled in an all-to-all architecture [1]. Although it is elegantly 
simple, the Kuramoto model is sufficiently flexible to be 
adapted to many different contexts, hence it is widely used 
and is regarded as one the most representative models of cou- 
pled phase oscillators [2]. Recently, the Kuramoto model has 
received increased attention. For example, the authors in [3], 
[4], [5] discussed synchronization conditions for the Kuramoto 
model. The work in [6] gave a synchronization condition for 
delayed Kuramoto oscillators. Results are also obtained for 
Kuramoto oscillators with coupling topologies different from 
the original all-to-all structure. For example, the authors in [7] 
and [8] considered the phase locking of Kuramoto oscillators 
coupled in a ring and a chain, respectively. Using graph theory, 
the authors in [9], [10], [11] discussed the synchronization of 
Kuramoto oscillators with arbitrary coupling topologies. The 
authors in [12] proved that exponential synchronization can be 
achieved for Kuramoto oscillators when phases lie in an open 
half-circle. 

Studying the influence of the pacemaker (also called the 
leader, or the pinner [13]) on Kuramoto oscillators is not 
only of theoretical interest, but also of practical importance 
[14], [15]. For example, in circadian systems, thousands of 
clock cells in the brain are entrained to the light-dark cycle 
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[16]. In the clock synchronization of wireless networks, time 
references in individual nodes are synchronized by means of 
intercellular interplay and external coordination from a time 
base such as GPS [17]. Hence, Kuramoto oscillators with a 
pacemaker are attracting increased attention. The authors in 
[14] and [18] studied the bifurcation diagram and the steady 
macroscopic rotation of Kuramoto oscillators forced by a pace- 
maker that acts on every node. Based on numerical methods, 
the authors in [19] showed that the network depth (defined 
as the mean distance of nodes from the pacemaker, a term 
closely related to pinning-controllability in pinning control 
[20]) affects the entrainment of randomly coupled Kuramoto 
oscillators to a pacemaker. Using numerical methods, the au- 
thors in [21] discovered that there may be situations in which 
the population field potential is entrained to the pacemaker 
while individual oscillators are phase desynchronized. But 
compared with the rich results on pacemaker-free Kuramoto 
oscillators, analytical results are relative sparse for Kuromoto 
oscillators forced by a pacemaker. And to our knowledge, there 
are no existing results on the synchronization rate of arbitrarily 
coupled Kuramoto oscillators in the presence of a pacemaker. 

The synchronization rate is crucial in many synchro- 
nized processes. For example, in the main olfactory system, 
stimulus-specific ensembles of neurons synchronize their firing 
to facilitate odor discrimination, and the synchronization time 
determines the speed of olfactory discrimination [22]. In 
the clock synchronization of wireless sensor networks, the 
synchronization rate is a determinant of energy consumption, 
which is vital for cheap sensors [23], [24]. 

We consider the exponential synchronization rate of Ku- 
ramoto oscillators with an arbitrary topology in the presence 
of a pacemaker. In the identical natural frequency case, we 
prove that synchronization (oscillations with identical phases) 
can be ensured, even when phases are not constrained in 
an open half-circle. In the non-identical natural frequency 
case where perfect synchronization has been shown cannot be 
achieved [2], [25], we prove that phase locking (oscillations 
with identical oscillating frequencies) can be ensured and 
synchronization can be achieved in the sense that phase 
differences can be reduced to an arbitrary level. In both cases, 
the influences of the pacemaker and local coupling strength 
on the synchronization rate are analyzed. 

II. Problem formulation and Model 

TRANSFORMATION 

Consider a network of N oscillators, which will hence- 
forth be referred to as 'nodes'. All N nodes (or a subset) 
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receive alignment information from a pacemaker (also called 
the leader, or the pinner [13]). Denoting the phases of the 
pacemaker and node i as (po and fi, respectively, the dynamics 
of the Kuramoto oscillator network can be written as 



ifO = Wq 



N 



(pi = W t + ^ a t,j sin (<Pj ~ fi) + 9i sin (<Po - fi) 



(1) 



for 1 < i < N, where u>o and Wi are the natural fre- 
quencies of the pacemaker and the ith oscillator, respectively, 
ctij sin (ifj — (fi) is the interplay between node i and node j 
with Q>i j > denoting the strength, gi sin(y?o — fi) denotes 
the force of the pacemaker with g L > denoting its strength. 
If dij = (or gi = 0), then oscillator i is not influenced by 
oscillator j (or the pacemaker). 

Assumption 1: We assume symmetric coupling between 
pairs of oscillators, i.e., djj = aji. 

Next, we study the influences of the pacemaker, gi, and local 
coupling, ciij, on the rate of exponential synchronization. 

Solving the first equation in ([TJ gives the dynamics of the 
pacemaker fo = w^t + <fro, where the constant </>o denotes 
the initial phase. To study if oscillator i is synchronized to 
the pacemaker, it is convenient to study the phase deviation of 
oscillator i from the pacemaker. So we introduce the following 
change of variables: 



fi 



fa + Ci = w ot + <po + Ci 



(2) 



Ci E [— 2tt, 2tt] denotes the phase deviation of the ith oscillator 
from the pacemaker. Due to the 27r-periodicity of the sine- 
function, we can restrict our attention to Ci € [— ir, ir}. 
Substituting (O into (fl} yields the dynamics of Ci- 



Ci 



Wi 



w 



N 



; sin (Cj - Ci) ~ 9i sin(&) (3) 



Since Ci is the relative phase of the ith oscillator with respect 
to the phase of the pacemaker, it will be referred to as relative 
phase in the remainder of the paper. 

By studying the properties of (01, we can obtain: 

• Condition for synchronization: If all Ci converge to 0, 
then we have fx = f 2 = ■ ■ ■ = fN = fa when t — > oo, 
meaning that all nodes are synchronized to the pacemaker. 

> Exponential synchronization rate: The rate of synchro- 
nization is determined by the rate at which Ci decays to 
0, namely, it can be measured by the maximal a satisfying 

\\C(t)\\<Ce- at \\C(0)l C = [Ci, C2, ■ ■ ■ , Cn] T (4) 

for some constant C, where j| • j is the Euclidean norm. 
a measures the exponential synchronization rate of (O: 
a larger a leads to a faster synchronization rate. 
Remark 1: When Wi and wq are non-identical, synchroniza- 
tion (Ci = 0) cannot be achieved in general. But we will 
prove in Sec. HV-Cl that the synchronization error can be made 
arbitrarily small by tuning the strength of the pacemaker gi. 

Assigning arbitrary orientation to each interaction, we can 
get the N x M incidence matrix B (M is the number of 
interaction edges, i.e., non-zero ajj (1 < i < N, j < i)) of 



the interaction graph [26]: _B,; j — 1 if edge j enters node i, 
-1 if edge j leaves node i, and B t j = otherwise. 



B 



Then using graph theory, (O can be recast in a matrix form: 



C = - Gsin£ - BWsm. (B T C) 



(5) 



where = [wi — w , W2 — w , . . . , w^ — Wo] T , G = 
diag(gi, g 2 , ffjv)> and W = diag(^i, v> 2 , . . . , v M ). Here 
i/j (1 < i < M) are a permutation of non-zero cii j (1 < i < 
N,j < i) and diag(») denotes a diagonal matrix. 



III. The identical natural frequency case 
When w\ = w% = . . . = wn = Wo, © reduces to: 
C= -G sin£ - BW sin (B T £) 



(6) 



To study the exponential synchronization rate, we first give 
a synchronization condition: 

Theorem 1: For the network in d6]l, denote e = max \d\ 

l<i<N 

and sinc(x) = sin(x)/a;, then 

1) when £ < §, the network synchronizes if at least one gi 
is positive and the coupling ajj is connected, i.e., there 
is a multi-hop link from each node to every other node; 

2) when ^ < e < ir, the network synchronizes if the 
following inequality is satisfied: 



> 



max 



J siiic(2£ )A max (£W.B T ) 
— sinc(£) 



N 

1 ^ , sin(e) 



(7) 



where A max (») denotes the maximal eigenvalue, g m - ln 
and £o G , 7r) are determined by 

5min = min{5i, . . . , g N }, 2s a cos (2e ) = sin(2£ ) (8) 
Proof: We first prove that when C £ [ — e, e] x . . . x 
[— £, e] = [— £, e] N where x denotes Cartesian product, they 
will remain in the interval under conditions in TheoremQ] i.e, 
the n-tuple set [—e, e] N is positively invariant for ©. 

To prove the positive invariance of [—£, e] N , we only need 
to check the direction of vector field on the boundaries. When 
£ < |, if Ci = £> we have —ir < — 2e < Cj ~ Ci < f° r 
1 < j < N. So in ®, sin(^ - &) < Q and sin(^) > 
hold, and hence Ci < holds (Note that Wi — w$ = 0). Hence 
the vector field is pointing inward in the set, and no trajectory 
can escape to values larger than e. Similarly, we can prove that 
when Ci = Ci > holds. Thus no trajectory can escape 



to values smaller than —e. Therefore 



-£, £ 



N 



is positively 



invariant when £ < f ■ When j < e < ir, if Ci = s, we have 
sin£i = sin£ > and sin(£j — Ci) < 1 f° r 1 < j < N. 
So when Wi — wq, if © is satisfied, the right hand side of 
(O is negative, i.e., < holds. Therefore the vector field 
is pointing inward in the set and no trajectory can escape 
to values larger than e. Similarly, we can prove that if Ci = 
— £, Ci > holds under condition (0. Thus no trajectory can 
escape to values smaller than — e. Therefore [—e, e] is also 
positively invariant for ? < e < tt if is satisfied. 

Next we proceed to prove synchronization. Define a Lya- 
punov function asV = ±C T C- V >0 is zero iff all Ci are zero, 
meaning the synchronization of all nodes to the pacemaker. 
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(9) 



(10) 



Differentiating V along the trajectories of (|6]l yields 

V = £, T £, = -£ T (G sin £ + BW sin(B T £)) 
= -£ T GSi£ - £ T BWS 2 B T £, 

where Si € TZ NxN and S 2 G K AIxM are given by 

51 = diag {sinc(£i), . . . , sinc(£ w )} , 

5 2 = diag {sinc(S T 0i, • • • , sinc(.B T £)M} 

with (B T £)i denoting the zth element of the Mxl dimensional 
vector B T £. 

From dynamic systems theory, if GSi + BWS 2 B T in (O is 
positive definite when £ 7^ 0, then F is negative when £ 7^ 
and 1/ will decay to zero, meaning that £ will decay to zero 
and all nodes are synchronized to the pacemaker. 

1) When all & are within [-e, e] with < s < § , (B T £) 4 is 
in the form of £ m — £„ (1 < m, n < iV), and hence is restricted 
to (— 7T, tt). Given that in (— tt, tt), sinc(x) > holds, it 
follows that Si and S2 satisfy the following inequalities: 



5*i > ciJ, (Ji = min sinc(ai) = sinc(e), 

— e<.T<e 

S*2 > 02-f, 02 — min sinc(a;) = sinc(2e) 

-2e<x<2e 



(ID 



So we have GSi + BWS 2 B T > a x G + a 2 BWB T , which in 
combination with (O produces 

V < -£ T (aiG + a 2 BWB T ) £ (12) 

It can be verified that <7iG + a 2 BWB T is of form: 

(TiG + a 2 BWB T = (Tidiag{gi, g 2 , . . . , g N } + o 2 L (13) 

with L € ji NxN given as follows: for i 7^ j, its (z, j)th element 



-dij, for i = j, its (i,j')th element is 



JV 



Since cti and a 2 are positive, and a^j are non-negative, 
it follows from the Gershgorin Circle Theorem that <7iG + 
a 2 BWB T only has non-negative eigenvalues [27]. Next we 
prove its positive definiteness by excluding as an eigenvalue. 

Since the topology of Ojj is connected, a\G + a 2 BWB T 
is irreducible from graph theory [27]. This in combination 
with the assumption of at least one gi > guarantees that 
<JiG + (T 2 BWB T is irreducibly diagonally dominant. So from 
Corollary 6.2.27 of [27], we know its determinant is non-zero, 
and hence is not its eigenvalue. Therefore <JiG + a 2 BWB T 
is positive definite, and hence V will converge to 0, meaning 
that the nodes will synchronize to the pacemaker. 

2) When & G [-e, e] (1 < i < N) with f < e < tt, 
Si is positive definite but S2 is not since {B T £)i is in 
[— 2e, 2e], and thus sinc(£? T £)i may be negative. It can be 
proven that sinc(x) is monotonically decreasing on [0, 2eo] 
and monotonically increasing on [2eo, 27r] (using the first 
derivative test), where £0 £ (f , tt) is determined by ([8j. 
Hence we have Si > sinc(e)/ and S2 > sinc(2e )I where 
sinc(2e ) < 0. Therefore (O reduces to 



Thus £ 
holds. 



V < -sinc(£)£ T G£ - sinc(2e )£ T BWB T i 
< -£ T (sinc(e)G + sinc(2£ )BI¥ J B T ) £ 

-> if 5 min sinc(e) + sinc(2e )A max (BVFS T ; 



(14) 



> 



Remark 2: It is already known that for general Kuramoto 
oscillators without a pacemaker, synchronization can only be 
ensured when max cpi — min tpi is less than tt, i.e., the initial 
phases lie in an open half-circle [9], [10], [11], [12], [28], [29] 
(although when phases are lying outside a half-circle, almost 
global synchronization is possible by replacing the sinusoidal 
interaction function with elaborately designed periodic func- 
tions [30], [31], it may introduce numerous unstable equilibria 
[31]). Here, synchronization is ensured even when ^ is outside 
(— f , § ), i.e., when phase difference (fli — tpj = £»— £j is larger 
than tt, meaning that the phases can lie outside a half-circle. 
This shows the advantages of introducing a pacemaker. 

Remark 3: Theorem Q] indicates that when £j is outside 
(— f , § ), i.e., when phases cannot be constrained in one open 
half-circle, all nodes have to be connected to the pacemaker 
to ensure synchronization. In fact, when some oscillators are 
not connected to the pacemaker, the relative phases may 
not converge to 0. For example, consider two connected 
oscillators, 1 and 2, with coupling strength ai <2 = &2.1 = «• If 
the pacemaker only acts on oscillator 1 with strength gi = k 
and the phases of the pacemaker, oscillator 1 and 2 are tt, OAtt, 
and 1.67T, respectively, though £1 = — 0.6-7T and £2 = 0.67T are 
all within [— 0.6-7T, O.671-], numerical simulation shows that £2 
will not converge to no matter how large k is. 

Remark 4: Since the eigenvalues of BWB T are non- 
negative [26], X max (BWB T ) > 0. 

Based on a similar derivation, we can get a bound on the 
exponential synchronization rate: 

Theorem 2: For the network in d6]l, denote e = max 

l<i<JV 

If the conditions in Theorem[T]are satisfied, then the exponen- 
tial synchronization rate can be bounded as follows: 

1) when < e < § holds, the exponential synchronization 
rate is no worse than 

ai = min {£ T (a x G + a 2 BWB T ) £/(£ T 0} 

? (15) 
= A min (aiG + a 2 BWB T ) 

with ctiG + a 2 BWB T given in ( Tf3l : 

2) when f < e < 7r holds, the exponential synchronization 
rate is no worse than 

a 2 = 5 min sinc(£) + sinc(2£ )A max (SW / B T ) (16) 
Proof: From the proof in Theorem [T] when < 
e < f, we have V < —2aiV, which means V(t) < 
C 2 e- 2ait V(0) =S> ||£(t)|| < Ge- Qlt |j£(0)|| for some positive 
constant G. Thus the synchronization rate is no less than ai. 

Similarly, when | < e < tt holds, we have V < ~2a 2 V. 
Hence the exponential synchronization rate is no less than a 2 , 
which completes the proof. ■ 

Remark 5: When < e < \ holds and there is no 
pacemaker, i.e., G = 0, using the average phase (p = 
2~2i=i at as re f erence 7 we can define the relative phase as 
& = & - <p. Since £ T 1 = with 1 = [1, 1, . . . , 1] T , 
the constraint £ T 1 = is added to the optimization 
min{£ T (aiG + a 2 BWB T )C/(C T 0) in G3- Given that 

G = and BWB T is the Laplacian matrix of interaction 
graph and hence has eigenvector 1 with associated eigenvalue 
[27], A m in in (fT3T > reduces to the second smallest eigenvalue, 
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which is the same as the convergence rate in section IV of [32] 
obtained using contraction analysis. 

Eqn. (fT6] l shows that when max|^| = s > j, a stronger 

i 

pacemaker, i.e., a larger g m i n leads to a larger a 2 , but the 
relation between ol\ and gi when max|^| = e < ^ is not 
clear. (In this case, g mm may be zero.) We can prove that in 
this case a± also increases with gi for any £ = 1,2,..., JV: 

Theorem 3: Both ot\ in (fT3T > and «2 in ( fl~6T > increase with 
an increase in pacemaker strength. 

Proof: As analyzed in the paragraph above Theorem 
[3j we only need to prove Theorem [3] when e < ? holds, 
i.e., ai is an increasing function of Recall from $13[ 
that <j\G + a 2 BWB T is an irreducible matrix with non- 
positive off-diagonal elements, so there exists a positive p 
such that pi — (a\G + a 2 BWB T ) is an irreducible non- 
negative matrix. Therefore, A max (pi — (o~iG + a 2 BWB T )) 
is the Perron-Frobenius eigenvalue of pi — (<7iG + 
a 2 BWB T ) and is positive [27]. Given that for any 1 < 
i < N, fi — \i(a\G + a 2 BWB T ) is an eigenvalue 
of pi — (<7iG + a 2 BWB T ) where A^ denotes the ith 
eigenvalue, we have p — A min (cr 1 G + a 2 BWB T ) = 
A max (a«/ - (a 1 G + a 2 BWB T )), i.e., oi = A min (<7iG + 
a 2 BWB T ) = fi- A max (pi - (ctjG + a 2 BWB T )). 

Since the Perron-Frobenius eigenvalue of pi — (a\G + 
a 2 BWB T ) is an increasing function of its diagonal elements 
[27], which are decreasing functions of all it follows that 
Amax (l^I — (ciG + a 2 BW B T )} is a decreasing function of 
gi, meaning that ai is an increasing function of all gi. ■ 

Remark 6: When all & are in (— f, §), since S 2 in (fTOb 
is positive definite, which leads to —^ T BWS 2 B T ^ < 0, the 
local coupling will increase cti in <TT3T >. But when max|£j| is 

i 

larger than ^, S 2 can be indefinite, hence -£ T BWS 2 B T £, 
can be positive, negative or zero, thus the local coupling may 
increase, decrease or have no influence on the synchronization 
rate. This conclusion is confirmed by simulations in Sec. [V] 

IV. The non-identical natural frequency case 

When natural frequencies are non-identical, Kuramoto os- 
cillators cannot be fully synchronized [2], [25]. Next, we will 
prove that synchronization can be achieved in the sense that the 
synchronization error (defined as the maximal relative phase) 
can be made arbitrarily small. This is done in two steps: first 
we show that under some conditions, the oscillators can be 
phase-locked, then we prove that the relative phases can be 
trapped in [—6, 5} for an arbitrary <5 > if the pacemaker 
is strong enough. The role played by the phase trapping 
approach is twofold: on the one hand, it makes the conditions 
required in phase locking achievable, and on the other hand, 
in combination with the phase locking, it can reduce the phase 
synchronization error to an arbitrary level. 

A. Conditions for phase locking 

When the natural frequencies are non-identical, the dynam- 
ics of the oscillator network are given in (|5). As in previous 
studies, we assume that the natural frequencies are constant 
with respect to time. The results are summarized below: 



2) f 



Cli j COS 



(2e) 



Theorem 4: Denote e = max|£j|, then the network in (0 

i 

can achieve phase locking if 

1) < e < ? holds, at least one gi is positive, and the 
coupling cii j is connected; 

< £ < f and g min 

;os(2e)X^BWB-) r g 

cos e 7 • I "t— < 

hold. 

Proof: To prove phase locking, i.e., all oscillators oscil- 
late at the same frequency, we need to prove that the oscillating 
frequencies ipi are identical. From (0, we have 99, = wo + £i, 
so if C = £ converges to zero, then phase locking is achieved. 
Differentiating © yields 

C = -GS 3 C ~ BWSaB t Q (17) 

where 

53 = diag(cos£i, cos£ 2 , • ■ • , cos^v) , 

5 4 = diag (cos(B T £) 1 , cos(B T £) 2 , . . . , cos(B T £) M ) 

Following the line of reasoning of the proof of Theorem Q] 
we can prove that ( is positively invariant under conditions in 
Theorem [4] Next we proceed to prove the convergence of £. 

Define a Lyapunov function as V = |C T C- Differentiating 
V along the trajectory of ( fTTI i yields 

V = C T C = -C T GS 3 C - ( T BWS 4 B T <; (19) 

Following the line of reasoning of Theorem Q] when ( 7^ 0, 
we can obtain V < under the conditions in Theorem |4] So 
V, and hence ( will converge to 0. Thus oscillating frequencies 
become identical and phase locking is achieved. ■ 

Remark 7: In the absence of a pacemaker, the authors in [3] 
proved that if the phase difference between any two oscillators, 
i.e., tpi — ifj, Vi, j, is within [— f , § ], then phase locking can 
be achieved. Given tpi — tpj = Vi,j, the condition in [3] 

only applies to — ? < £j < ? in our formulation framework. 



B. A bound on the exponential rate of phase locking 

Theorem 5: For the network in ©, denote e = max If 

. . i 

the conditions in Theorem |4] are satisfied, then 

1) when < e < j holds, the exponential phase-locking 
rate is no worse than 



a 3 Ami„ {<t 3 G + a A BWB T ) 



(20) 



with (73 = cos e and 04 = cos 2e; 
2) when ? < e < ? holds, the exponential phase-locking 
rate is no worse than 

c*4 = ffmincos(e) + cos(2e) \ max (BWB T ) (21) 
Proof: Theorem [5] can be derived following the line of 
reasoning of Theorem [2] and thus is omitted. ■ 
Remark 8: Following Theorem [3] we can prove that a 
stronger pacemaker always increases 03 (and 04). But a 
stronger local coupling can have different impacts: when < 
e < f , S 4 in flS is positive definite, -( T BWS 4 B T ( is 
negative, so the local coupling will increase 03. However, 
when j < e < ?, since S4 in (fT~ST > can be indefinite, 
— £ T BWS4B T ( can be positive or negative. Thus the local 
coupling may increase or decrease the rate of phase locking. 
The conclusion will be confirmed by simulations in Sec. [V] 
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C. Method for trapping relative phases 

In this section, we will give a method such that the relative 
phases are trapped in any interval [—6, S] with an arbitrary 

< S < 7T. 

Theorem 6: For © with frequency differences Q, denote 
e = max and ||0|| = V^WTl, then the relative phases can 

i 

be trapped in a compact set [—6, S] for an arbitrary < 5 < tt 
1) if < e < f and the following condition is satisfied: 

> ||0||/(<5sinc(e)) (22) 



9n 

2) if ^ < e < 7r and the following condition is satisfied 

smc(2eo)A max (£WB T ) 



5min > ||ft||/(<Ssinc(e))- 



(23) 



sinc(e) 

where Eq is defined in (|8}. 
Proof: Differentiating Lyapunov function V = ^£, T £, 
along the trajectory of © yields 



v=et=t 



£ T G sin £ - sin(.B T £) 



(24) 



= £ T tt - eGS^ - eBWS 2 B T £, 

with Si and S% defined in ( fTob . 

1) When < £ < f holds, we have Si > sinc(e)/ > 



and 52 > from previous analysis. Using d24l i. ( [Tot , 
and the fact A m i n (G) = <? m i n , we have 

V^< ||e||||0||- gmin sinc( £ )||e|| 2 (25) 
If & is outside [—6, 5} for some i, we have ||£|| = 
\Jj2iLi Q > ^> which in combination with ( f22t leads 
to V" < 0. Therefore all & will converge to [—5, 5]. 
2) When f < £ < 7r holds, from the analysis in Theorem 
[TJ we have Si > sinc(e)/ > and S% > sinc(2e )7. 
Then using d24l ) and the fact A m i n (G) = .g m in, we have 

v< Hellion - 5min sinc( e )ii£ii 2 

- smc(2e )X m ^(BWB T )U\\ 2 



(26) 



If £j is outside [—6, 5} for some i, we have ||£| > 6, 
which in combination with d23l leads to F < 0. Thus 
all & will converge to the interval [—5, <5]. ■ 
Remark 9: Theorem |6] used the important fact that if ||£|| = 

^IJi is restricted to the interval [0, 5], then all 6 are 
restricted to the interval [—(5, J]. 

Remark 10: When ||^|| < 7r, [3] gives a condition under 
which can be trapped in an arbitrary compact set. Since for 

a large number of oscillators N, ||£|| = \Jj2iLi < tt is 
difficult to satisfy, our result is more general. 

V. Simulation results 

We consider a network composed of TV = 9 oscillators. The 
coupling strengths a,,j are randomly chosen from the interval 
[0, 0.1]. They were found to form a connected interaction 
graph. As in previous studies, we use the modulus of the 



order parameter r 



N 2-J%= 



to measure the degree 



of synchrony [25]. The value of r (r € [0, 1]) will approach 1 
as the network is perfectly synchronized, and if the phases 



are randomly distributed [25]. According to [25], we have 
r w 1 when the oscillators are synchronized. So we define 
synchronization to be achieved when r exceeds 0.99. 

When the natural frequencies are identical, we set the 
phase of the pacemaker <po to ipo = wot with Wq = 1 and 
simulated the network using initial phases tp i = <p + £ 4 
with G (— f , f) and initial phases tp i = ip + £ 4 with 
£i G (— 7r, tt), respectively. In the former case, we connected 
the first oscillator to the pacemaker and set gi = g, g 2 = 
93 = ... = 99 = 0. In the latter case, we connected all 
oscillators to the pacemaker and set gi = g% = ■ . . = gg = g. 
In both cases, we set g = 1. To show the influences of 
the pacemaker on the synchronization rate, we fixed a^j and 
simulated the network under different pacemaker strengths 
m x g, where m = 1,2,..., 10. To show the influences 
of local coupling on the synchronization rate, we fixed the 
strength of the pacemaker to 3g and simulated the network 
under different local coupling strengths m x aj j for all a iy j, 
where m = 1,2,..., 10. All the synchronization times are 
averaged over 100 runs with initial in each run randomly 
chosen from a uniform distribution on ( — ^, f ) (in the former 
case) or on (— tt, tt) (in the latter case). The results are given 
in Fig. 1. It is clear that a stronger pacemaker always increases 
the synchronization rate, whereas the local coupling increases 
the synchronization rate when all £j are within (— §, §), and 
it may increase or decrease the synchronization rate when the 
maximal/minimal £j is outside (— 
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Fig. 1 . Times to synchronization under different strengths of pacemaker/local 
coupling (with all oscillators having identical natural frequencies). 
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Fig. 2. Times to phase locking under different strengths of pacemaker/local 
coupling (with oscillators having non-identical natural frequencies). 

When the natural frequencies are non-identical, we simu- 
lated the network using initial phases (fi = tpo + £j with 
& e (-J > f ) and initial phases ^ = ip + & with & G 
(— f , § ), respectively. In the former case, we connected the 
first oscillator to the pacemaker and set gi = g. In the latter 
case, we connected all the oscillators to the pacemaker and 
set 5i = <?2 = • • ■ = .99 = 9- The natural frequencies were 
randomly chosen from (0, 1). Tuning the strengths in the same 



6 



way as in the identical natural frequency case, we simulated 
the network under different strengths of the pacemaker and 
local coupling. All of the times to phase locking are averaged 
over 100 runs with initial £j randomly chosen from a uniform 
distribution on {—\, f ) (in the former case) or on (— f , §) 
(in the latter case). The results are given in Fig. 2. It is clear 
that a stronger pacemaker always increases the rate to phase 
locking, whereas the local coupling increases the rate to phase 
locking when all & are within (— j , j), and it may increase or 
decrease the rate to phase locking when the maximal/minimal 
& is outside (-§ , §). 

To confirm the prediction that can be made smaller by 
making the pacemaker strength stronger, we set g% = . . . = 
gg = g and simulated the network under initial phases ipi = 
<Po + £,i with £i G (— 7r, 7r). Using the same £j, the maximal 
final relative phase when the strength of the pacemaker g is 
made to (to = 1,2, ...,10) times greater is recorded and 
given in Fig. [3] It can be seen that the maximal final relative 
phase (i.e., synchronization error) decreases with the strength 
of the pacemaker, confirming the prediction in Theorem [6] 
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0.45 




2 4 6 8 10 
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Fig. 3. The maximal final relative phase (phase synchronization error) 
under different strengths of the pacemaker when oscillators have non-identical 
natural frequencies (which are randomly chosen from the interval (0, 1)). 

VI. Conclusions 

The exponential synchronization rate of Kuramoto oscilla- 
tors is analyzed in the presence of a pacemaker. In the identical 
natural frequency case, we prove that synchronization to the 
pacemaker can be ensured even when the initial phases are 
not constrained in an open half-circle, which improves the 
existing results in the literature. Then we derive a lower bound 
on the exponential synchronization rate, which is proven an 
increasing function of the pacemaker strength, but may be 
an increasing or decreasing function of the local coupling 
strength. In the non-identical natural frequency case, a similar 
conclusion is obtained on phase locking. In this case, we also 
prove that relative phases (synchronization error) can be made 
arbitrarily small by making the pacemaker strength strong 
enough. The results are independent of oscillator numbers in 
the network and are confirmed by numerical simulations. 
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